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_^ ' We present a first calculation of the strange quark mass using domain wall fermions. This 
paper contains an overview of the domain wall discretization and a pedagogical presentation 

fvq ' of the perturbative calculation necessary for computing the mass renormalization. We combine 

f^ , the latter with numerical simulations to estimate the strange quark mass. Our final result in 

0^ ' the quenched approximation is 95(26) MeV in the MS scheme at a scale of 2 GeV. We find that 

^j , domain wall fermions have a small perturbative mass renormalization, similar to Wilson quarks, 

"*^ ' and exhibit good scaling behavior. 
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1 Introduction 

The determination of the quark masses from first principles is an important task facing particle 
theorists today. The light quark masses are among the most poorly determined parameters of the 
Standard Model. The cause of the difficulty is the confining nature of QCD: quarks exist only in 
bound states. Furthermore, most of the mass of light hadrons is due to the energy of the color fields 
surrounding the quarks rather than the quarks themselves; therefore a nonperturbative treatment 
of QCD is required to connect the quark masses in the QCD Lagrangian with the spectrum of 
hadronic states measured experimentally. 

Next-to-lowest order chiral perturbation theory (xPT) quite precisely predicts the ratios of 
quark masses [Q] but cannot set the absolute scale. The most promising method of computing the 
light quark masses (i.e., m^, m^ and mg) is lattice QCD. It has been suggested that QCD sum 
rules can be used to place fairly strict lower bounds on mg + tUu and rud + 7n„ by using analyticity 
conditions Q; however, calculations of the values of the light quark masses from sum rules are 
thought to involve many uncertainties ||3|. 

The feasibility of calculating the quark masses through Monte Carlo simulation of lattice QCD 
has been recognized since the early days of the field Q. Most previous work utilized two formula- 
tions of lattice fermions: Wilson fermions which explicitly break chiral symmetry at finite lattice 
spacing and suffer from large discretization errors; and Kogut-Susskind fermions which maintain 
a remnant chiral symmetry, but badly break flavor symmetry and seem to have poorly converged 
weak-coupling expansions. Recently, Sheikholeslami-Wohlert (SW) fermions ||5|, an improvement 
of Wilson fermions, have also been used to compute the light quark masses |6|, |^. 

The usual method of computing the light quark masses on the lattice is the following. For fixed 
gauge coupling and various bare quark masses, one computes the pseudoscalar or vector meson mass 
from the exponential decay in Euclidean time of an appropriate correlation function. Using leading 
order xPT one extrapolates in the bare quark mass niq to where the meson mass takes on its physical 
value, then the corresponding value of the bare quark mass can be converted to the renormalized 
quark mass in any desired scheme. This method implicitly uses the vector Ward identity, so the 
difference rUg — rric is renormalized by Z^. = Zg , where mc is the bare quark mass corresponding 
to zero pion mass and Zg is the renormalization constant of the scalar density. Usually one converts 
the renormalized lattice quark mass to the continuum MS regularization scheme by matching the 
weak coupling expansions of Zg in both schemes. Using this procedure the light quark mass 
mi = {jriu + m(i)/2 and the strange quark mass rus may be computed independently, for example, 
using the pseudoscalar spectrum to fix mi and the vector spectrum to fix m^. A comprehensive 
analysis of the light quark masses using this method appears in Ref. [Q. 

Recently several attempts have been made to remove some sources of uncertainty in the usual 
method. A large source of error in Wilson fermion calculations is the determination of the chiral 
limit; since they explicitly break chiral symmetry, Wilson quarks become massless at a nonzero 
critical bare quark mass, ?n,c 7^ 0. One can avoid this error by using the axial Ward identity to 
fix the bare quark mass ||9|, |l^]. One computes (0|9o^oK) ^ind (0|P|7r), where j4^ is the local non- 
singlet axial vector current and P the non-singlet pseudoscalar density, and then the quark mass 
is given by the ratio (0|5oAo|7r)/(0|P|7r) and is renormalized by Zm = Za/Zp. Since the vector 



meson and baryon spectra cannot be used with this method, only one of either mi or m^ may be 
fixed independently, the other is necessarily related by chiral perturbation theory. 

Another large uncertainty enters into the matching between lattice and continuum regulariza- 
tions. The typical scale at which this matching occurs is 2 GeV where the validity of WCPT is 
tenuous. In the case of Kogut-Susskind fermions, lattice WCPT is untrustworthy: next -to-leading 
order corrections can be 50-100% of the leading order term. Therefore, nonperturbative calculation 
of the renormalization factors is very desirable. Two methods are being explored presently which 
may remove the need for a perturbative expansion of the lattice theory lO], or push it to a very 
high energy scale at which the expansion parameter is much smaller |jl2|. Finally, the quenched 
approximation seems to give quark masses which are roughly 20% larger than unquenched quark 
masses [^]. Clearly this indicates that full QCD simulations are necessary for a precise calculation 
of light quark masses. 

Presently no consensus has been reached regarding the values of the light quark masses, even 
within the quenched approximation. For example, using Wilson fermions and perturbative match- 
ing, the strange quark mass is 115(2) MeV when the kaon is used to fix the bare quark mass and 
143(6) MeV when the (f) meson is used [|l3|, where the mass is defined in the MS scheme at 2 GeV.|^ 
Results using the SW action give a lighter strange quark mass of 95(16) MeV |0]. Furthermore, 
an exploratory nonperturbative determination of the quark mass renormalization agrees with the 
perturbative renormalization for the usual quark mass definition, but differs with the perturbative 



renormalization for the axial Ward identity definition [11]. Using the nonperturbative renormal 



ization and the axial Ward identity, Ref. ||ll[ finds a strange quark mass of 130(18) MeV. A more 



comprehensive presentation of the current status appears in Ref. |14]. 

In this paper, we employ a new fermion discretization to compute the light quark masses: do- 
main wall fermions. Domain wall fermions utilize a fictitious extra (in this case, fifth) dimension in 
order to preserve chiral symmetry at nonzero lattice spacing; the chiral symmetries of the contin- 



uum are exactly preserved in the limit of an infinite fifth dimension |18, 56]. The idea originated 



in the context of chiral gauge theories. In Ref. ]15], Kaplan constructed free lattice chiral fermions. 



without doublers, in 2k dimensions by considering Dirac fermions in 2k + 1-dimensions coupled 
to a mass defect in the extra dimension, or domain wall. For periodic boundary conditions in the 
extra dimension, an anti-domain wall also appears which supports 2A;-dimensional chiral fermions 
of the opposite handedness. Although the suitability of this approach for chiral gauge theories is 
still under intensive study, its usefulness for simulations of chirally symmetric vector gauge theories 
such as QCD now appears well established (for a review, see Ref. ]|lq] ). 

Since the first suggestion that domain wall fermions offer a way to study chiral symmetry 
breaking of QCD [0, |l^, considerable work has been done to assess the practicality of the method. 



In Ref. 117], a simplification of Kaplan's original proposal is made for QCD simulations: half of the 



extra dimension is discarded, and the domain walls effectively become the boundaries of the extra 



^ The widespread belief is that the quenched approximation yields a K — (j> splitting which is smaller than the 
experimental value. Using a regularization independent renormahzation scheme, recent quenched simulations with 
Kogut-Susskind fermions give m^ = 106(7) MeV with the kaon as input versus 129(12) MeV with the (j) as input [p4| . 
On the other hand, a recent quenched study using the SW action with a nonperturbatively determined coefficient 
claims to reproduce the physical K — (f> splitting if their chiral fit is quadratic rather than linear [pSl . 



dimension. For free field theory it has been shown that, for a range of the input parameters, a hght 
four-dimensional mode of definite chirality is bound to one boundary and a similar mode of opposite 
chirality is bound to the other boundary; the mixing between the two modes is exponentially 
suppressed with the size of the extra dimension |18, |^, |2^, 21]. Significant suppression of the 



mixing between the modes was also seen in nonperturbative simulations, but whether or not it is 



purely exponential remains an open question [g^, ^, ^, 16 1. Furthermore, the non-singlet axial 



Ward identity is reproduced, and predictions from chiral perturbation theory for the dependence of 



the pseudoscalar meson mass on the quark mass and for the kaon mixing parameter are satisfied |22, 



23, |l^. Also, the expected behavior of (qq) in the quenched approximation due to topological zero 



modes is reproduced |24, 25 1. 



The paper is structured as follows: Section y introduces the details of the domain wall fermion 
action. Section ^ contains the results of the one-loop calculation of the massive quark self-energy, 
and Section ^ gives the details of our Monte Carlo simulations. In Section |5| we combine analytical 
and numerical results to give a value for the strange quark mass. Finally, we present our conclusions 
in Section ^ and include some details of our calculation in the Appendix. 

2 Domain wall fermions 

In this Section we review some properties of domain wall fermions in QCD, following the original 
boundary fermion variant by Shamir [1^]. We take a pedagogical point of view and introduce 
notation and methods which will be relevant to our work later in this paper. We write down 
the action and the propagator, discuss the physics described by the light modes coupled to the 
boundaries, and mention previous work supporting the domain wall formulation of lattice QCD. 
For further details one should refer to the literature cited throughout the Section. 

2.1 The action 

On a lattice with spacing a, the domain wall fermion action is given by 

~"'^^Y1 i's{x/a)Ds^s' (x/a, y/a)ips' (y/a), (1) 

x,y s,s' 

where x, y are four-dimensional Euclidean spacetime coordinates and s, s' € [1, Ng] are coordinates 
in the fifth dimension. The Dirac operator can be separated into a four-dimensional part, D", and 
a one-dimensional part, D : 

aDs,s'{x,y) = aD^x,y) X 6s,s' + 6{x - y) x aDJ-.,,; (2) 

X and y have been rescaled to be dimensionless. The first term is the four-dimensional Wilson-Dirac 
operator with a mass term which is negative relative to the usual 4d Wilson fermion action: 



aDHx,y) = -J2[{l+^^)U^{x)6{x + il-y) + {l--f^)Ul{y)6{x-i:i-y) 

+ {aM-A)6{x-y). (3) 



For 1 < s < Ng, the Dirac operator in the extra dimension is given by 

aD^ „; = —6 



+ 2 (^ + 75)5s+i,s' + 2 (^ ~ 75)'^s-i,s'- 



(4) 



Note that the five-dimensional fermions are coupled to four-dimensional gauge fields which are 
identical at each s; i.e. the link matrices obey 



Uf,,s{x) = U^{x) for /ue [1,4], 
U^,s{x) = 1 for 1 < s < Ns. 



(5) 



The boundary conditions in the fifth dimension are anti-periodic with a weight am which, as we 
will see, is proportional to the Ad quark mass. The Dirac operator for the fifth dimension can be 
separated into its chiral components by the projectors, P± = (1 ± 75)/2, such that 



Di 



Ks'P+ + ^.:;^^- 



(6) 



In matrix notation, 



aD: 







1 
-1 1 





y —am 



\ 




-1/ 



aD 



-L,- 



/ -1 

1 -1 

1 

V 







1 



-am, \ 



-1 J 



(7) 



Since the perturbative calculation is simpler in (four-) momentum space, we Fourier transform 
the ordinary spacetime coordinates. The five-dimensional domain wall Dirac operator (|2|) becomes 



aDs^s'{ap) = [ ^ i7/^ sin ap^ + aM- ^(1 - cos ap^] 



+ aDi 



(8) 
(9) 



where ap^ = sinap^, and the W are related to the aD ' by 



W^ 



aDi'^ + 



aM — 2J(1 ~ cos ap^) 



(10) 



2.2 The mass matrix 

The discrete extra dimension can be interpreted as a flavor space with A^^ — 1 heavy fermions and 
one light flavor |26]. In that framework, the W in Eqn. (^) are mass matrices which govern the 
flavor mixing. One can have several physical pictures of how the mass hierarchy is maintained. The 
chiral symmetry is manifest in the domain wall picture; since the left- and right-handed modes 
are bound to opposite walls in the extra dimension and are separated by a distance Ng, so the 
chiral components can be rotated independently. The flavor space picture also provides insight. 



Ref. |27| relates the mass hierarchy in terms of a generaUzed see-saw formula. In fact, the Froggatt- 
Nielsen ^^ mechanism allows one to establish an approximate conservation law which protects the 
light mass from large radiative corrections |27| |. 

Let us examine the eigenvalues and eigenvectors of the tree-level mass matrix (in flavor space) 
for the action described above in Eqns. (|l|)-(^- Details are presented in Refs. ||l^, ^], and we 
repeat them in Appendix AJ so that they may be extended to the one-loop case. Since the mass 
matrix is not hermitian, we diagonalize the mass matrix squared. Let Q^ be the zero momentum 
limit of W- (0): 

/ —bo ... —am \ 
1 -6o ••• 
1 ... 



f^Os 



V 



(11) 



-bo J 



where 6o = 1 — M. Here and in the rest of the paper, we rescale M so that it is dimensionless. As 
shown in Refs. [17, ^, 19], when |6o| < 1 the smallest eigenvalue of riof^J (and of ilgOo) is 



(AW)2 = (am)2M2(2-M)2 + 0((am)^)+0((l-M)^= 



(12) 



Therefore, the mass of the light mode, given by A'^^, has an additive renormalization which is 
suppressed as A^^ — > oo for the range < M < 2. (We elaborate on this restriction on M in the 
next Section and throughout the remainder of the paper.) For |6o| < 1 the eigenvectors of ilor^Q 
are given by 

f JM{2 - M) e-"(^-i) (sign hoX~^ i = 1 
# = { ^ r r ,^r i\ (13) 

where the i— th eigenvector corresponds to eigenvalue (A*^*')^. The constant ao is defined through 

coshao = ^-^^- (14) 

Note that the light eigenmode in Eqn. ( |T^ ) is exponentially concentrated at the s = 1 boundary 
while the heavy (i > 1) modes are not. 

Once the eigenvectors of the mass matrix squared have been found, the Dirac operator can be 
diagonalized easily. Following Ref. pH] let us define unitary matrices U^^' and V^^^ such that 



Then the basis 



U. 

i^f-'Hp) 



(0) 



d') 



and yi°j ^ 



(0), 






(15) 






(16) 



diagonalizes DD^ and D^D. Furthermore, the Dirac operator itself is diagonal in this basis, up to 
terms which vanish as Nsexp{—aoNs) |^, ^ : 

= V^'^'^s(-p)(ia| + y(o)H^+C/(o)tp+ + C/WTy-y(o)tp_)v''^i^s(p); (17) 

where the y(o)pF+C/(°)t and C/(o)l^-y(o)t are diagonal in s,s' |0|. Let us define x = '4'f=i, the 
eigenstate of the lightest eigenvalue of the mass matrix (squared). This mode has the effective 
tree-level action 



ctree A 



■K/a ^Ap 



Since 



y(0)|y+[7{0)t 



=l,u=l 



/a (27r)4 



X(-P)(¥ + mM(2-M))x(p), 



(18) 



1,1 



a(^) m- 



2.3 The propagator 

The calculation of the tree-level fermion propagator Sf proceeds similarly to the diagonalization 
of the mass matrix presented above. The final expression for the propagator is complicated, so 



we write it here schematically; it is written explicitly for the present action in |17, |2^ and in 
Appendix BT. As in Ref. ^^ let us write the propagator as Sp = Dy{DD'^) and project out its 
chiral eigenstates so that 



sL'ip) 



where p^ = sinp^, and G (G ) is the inverse of D^D [DD^): 



G 



R 



p2 + w+W- 



(19) 



(20) 



We give G^ and G^ explicitly in Appendix B.l, but let us mention their general behavior. The 
homogeneous solutions of 



D^D 



Gs",s' — ^s 



and [DD 



(jr t.// a' 



(21) 



are exponentials: exp(zba(p)s). If < b{p) < 1 then the solutions are decaying exponentials; a{p) 
is real and defined through 

p2 + l + 62(p) 



where 



cosh a (p) 

Kp) 



2\b{p)\ ' 
l-M + ^(l-cosap^). 



(22) 
(23) 



If —1 < b{p) < then the solutions oscillate: (— 1)* exp(zba(p)s), with a as in Eqn. (p^). For 
\b{p)\ > 1 there is no longer a mode bound to the domain wall: the solutions of Eqn. ( pi| ) go like 
exp(zizia{p)s) where a is now defined through 

p^ + l + b\p) 
cos a{p) = ^^^ , (24) 

Thus, M must be in the range 

< M < 2 (25) 

in order for there to be a single massless fermion as A^^ ^ oo. At greater values of M, "doubler" 
states in the other corners of the Brillouin zone become nearly massless. For example, the four 
states with one component of momentum near n/a contribute, and so on in increments of 2 up to 
8 < M < 10 where again only a single state with all four momenta near n/a exists. In fact the 

action is symmetric under the changes M — > 10 — M and ips{x) — > (— l)^f ^'^ '^ipsix) so the physics 
of any region and the one reflected about M = 5 are identical ||2^. Henceforth, we concentrate 
our discussion on the single flavor theory near the origin of the Brillouin zone. Also note that the 
range of M in Eqn. ( p5|) is additively renormalized in the interacting theory, and we discuss this 
renormalization in detail in Sec. IST 



2.4 Summary of tree level properties 

Let us emphasize the following points. For strictly infinite Ng (and m = 0) there is a massless right- 
handed fermion bound to the s = 1 wall and a left-handed fermion bound to the anti-domain wall 



at s = oo |15]. For large, but finite, Ng the mixing between the two chiralities is exponentially 



suppressed at tree level |17]. One describes a light 4-component Dirac fermion by coupling the two 



light modes with an explicit chiral symmetry breaking term m [18|. Then the light Dirac fermion 



has a mass in the free theory equal to |17, 19 1 



am)^> = M{2 - M) am + {1 - My^' . (26) 

Therefore, neglecting exponentially small terms, domain wall fermions describe a light mode whose 
mass is multiplicatively renormalized. It also turns out that the light mode satisfies continuum-like 
axial Ward identities [^] . These features make domain wall fermions very attractive for simulating 
light quark physics, where chiral symmetry is crucial. 

Another virtue of the domain wall formulation is that in the limit Ng ^ oo the leading dis- 
cretization errors for the Ad effective action are 0{a'^). In the massless theory, any gauge invariant 
dimension-five operator with the required lattice symmetries can be written as a linear combination 
of the following two operators ||5|: 

Oi = qD^q (27) 

O2 = '^q<y^.uF^.q (28) 



where D^ is the second-order covariant derivative and F^^y = [D^^^Dy] is the field strength tensor. 
However, neither of these terms is invariant under a chiral transformation 

q -^ qe''^\ (29) 

Since chirahty violating effects have been shown to vanish as A^s — > oo, the contributions of Oi 
and O2 to the effective action must be suppressed. In this sense, the domain wall fermion action is 



an 0(a)-improved action |30, ^]. Even with an 0(a)-improved action, 0{a) errors can enter into 
observables as the operators may require 0{a) improvement as well. However, such improvements 
would also violate chiral symmetry and, by the same argument as above, are suppressed as Ns ^^ 00. 
Of course precise scaling tests are necessary to evaluate the extent of the improvement; although 



simulations to date [23|, including the ones in this work, are consistent with these expectations. 



3 Perturbative mass renormalization 

In this Section we present our one-loop calculation of the quark mass renormalization. A renor- 
malization factor is needed in order to match a lattice definition of quark mass to a continuum 
definition. We are also pursuing methods of computing renormalization factors nonperturbatively; 
however, that work is beyond the scope of this paper. 

We compute the matching factor Zm(/i,a) between a quark mass defined on a lattice with 
spacing a to a continuum quark mass renormalized at momentum scale /x: 

m^{^i) = Z„(/x,a)mLAT(a). (30) 

Zm is computed by equating the one-loop continuum fermion propagator to the one-loop lattice 
fermion propagator. In this Section we present the calculation of the full five-dimensional self- 
energy, and then we discuss its effect on M and m. The wavefunction renormalization has been 
computed already in Ref . [^0| ; we have extended that work to the massive case and present the full 
one-loop calculation here for clarity.|2| Only the main points are made in the body of the Section, 
while more details are given in Appendix ^. 

3.1 Five— dimensional fermion self— energy 

The fermion self-energy, S(p, ?7z), is given to one-loop order by the Feynman diagrams shown in 
Figure ||. We use p to denote the external momentum and I to denote the momentum in the loop 
integral. The tadpole graph has no fermion propagator in the loop, so it has trivial dependence on 
the fifth dimension, i.e. it is diagonal. On the other hand, the fermion in the loop of the half-circle 
graph may propagate in the fifth dimension (change flavor) while the gluon is unaffected. Therefore 
the half-circle graph has off-diagonal contributions in s,s' space. 

^ While this manuscript was in preparation, an independent one-loop calculation of the quark mass renormalization 
appeared in Ref. ISlI] which helped us in tracking down an error in our preliminary work 






(a) 






<i^ p-l '^a 



I 
(b) 






Figure 1: One-loop ferinion self-energy diagrams: the (a) tadpole and (b) half-circle graphs. 



Even with the extra dimension, the steps of evaluating the half-circle graph are much like those 
for the calculation using Wilson fermions [^, ^]. First an integral which has the same infrared 
(p -^ 0) limit is subtracted from S(p, m) to cancel logarithmic divergences. The difference may 
then be Taylor expanded about zero lattice spacing. In the continuum limit, we neglect terms in the 
expansion which vanish as a — > 0; since the coefficients of O^a) terms are exponentially suppressed 
with increasing A'^s the leading discretization errors are, in effect, 0{a}\ The resulting expression 
can be arranged as follows: 



-So + i^Si + mS- 



167r2 

After a lengthy calculation which is presented in Appendix [ 
into 



(31) 



these terms can be further subdivided 



y{0) _ 


_ r(0) 


- 2T6s 


v{i) 




+ c 




_ r(2) 


+ 43 



(32) 

T 

i^^s,s' (33) 

(34) 

where the L terms are proportional to In a, and the I and T terms are finite integrals to be computed 
numerically. We only need the renormalization of the lightest mode, so we delay further evaluation 
until we rotate the terms to the basis which diagonalizes the one-loop mass matrix. We should note, 
however, that the I and L terms are functions of M. Since M becomes additively renormalized, this 
dependence is the source of a systematic uncertainty: what numerical value of M should one use to 
compute I^^' and l'^)? We will address this issue in Sec. 3^. To summarize, the five-dimensional 
one-loop effective action is given by 



V's(-P) 



,-1 



+ 



W+P+ + w-p^ 
g^Cp 



+ 



1 + 



167r2 



9^Cf 
167r2 
T 
'2 



(j(0) 



2T 



(l« +/«-!)] +m|^(L(2)+/(2))}^^^^.,(p). (35) 



3.2 Renormalization of M 

We noted above (|25| ) that in the free theory M = is the critical point where light modes begin 
to appear bound to the domain walls. The 1/a contribution from the self-energy graphs, Sq given 



10 



in Eqn. (p2|), shifts this value in the same way that the niassless Hmit of Wilson-Uke ferniions is 
shifted: 

M = -^ M = Mc, (36) 

where Mc is the point at which the domain wall action first supports light chiral modes on the 
boundaries. If the tadpole contribution is dominant the shifts are identical since that graph is the 
same for domain wall and Wilson fermions. 

Of course, the whole range of M corresponding to light modes is renormalized when the coupling 
is non-zero, except that M = 5 is a fixed point since the action is still symmetric under M — > 10— M 
and ips{x) — > (— 1)"^^ ^'' ^ips{x). For example, for numerical simulations it is helpful to know the 
optimal value of M which minimizes the extent of the light mode in the fifth dimension. At 
tree-level the wavefunction corresponding to the zero mode bound to the s = 1 domain wall is a 
(5-function in the s-direction when M = 1; i.e., M = 1 is the optimal value for free domain wall 
fermions since the intrinsic quark mass arising from mixing of modes on opposite walls is minimized. 
However, simulations at /? ~ 6.0 have shown that the optimal M at which axial symmetries are 
preserved is somewhere around 1.7 and that axial symmetries are poorly respected at M = 1.0 

The question now becomes, how is the whole range of M renormalized? It cannot be a simple 
uniform shift since M = 5 is a fixed point. If we consider only the tadpole contribution to the 
self-energy, the shift is approximately uniform in each region (0 < M < 2, 2 < M < 4, . . .), 
weighted by a factor 

^(l-cos(p^))^ 4,2,0,-2,-4 (37) 

coming from the Wilson term. Surprisingly, this simple picture also describes the nonperturbative 
data well, as will be shown later in this paper. This tadpole-improved estimate of M was originally 
proposed in Ref. ||20|| . 

Perturbatively, the shift of Mc is given through the a~^ terms in Eqn. (|35|): 



P-^' = Ig li is not diagonal in the extra dimension, the one-loop calculation of Mc for the light mode 
involves rotating to the basis which diagonalizes the one-loop mass matrix.^ However, a reasonable 
first estimate for Mc assumes the tadpole graph is numerically much larger than the contribution 
from/(0): 

, = ^^, (39) 

where T is given in Eqn. ( ]56| ) and is numerically equal to 24.4. In Table IT] we give two values of 
M"^^, each computed with different definitions of the strong coupling constant, gy{?>Al/a) and 
gy{l/a). There is an obvious problem in deciding the relevant scale of this effect. In Section ^ 
we discuss the choice of coupling constant and scale in detail, but it is clear that a perturbative 

■^The calculation has been carried through in Ref. [ pl[ . 

11 



Table 1: Values of Merit estimated from the tadpole graph and computed with the 4-d Wilson 
fermion action. 





/3 = 6.3 


^ = 6.0 


(3 = 5.85 


Mt^'i(3.41/a) 
Mj:^'^{l/a) 


0.676 
0.985 


0.754 
1.17 


0.812 
1.41 


M^ 


0.1519 
0.708 


0.1572 
0.819 


0.1617 
0.908 



estimate of Mc is not precise enough for our purposes. We note, however, that for a reasonable 
estimate of the renormalized coupling, the shift due to the tadpole agrees well with the above 
nonperturbative estimate. 

Due to the above considerations, it is desirable to have a nonperturbative determination of 
Mc since it appears in the definition of the lattice quark mass. A direct search for Mc would be 
numerically prohibitive, so it is fortunate that there is a simple nonperturbative estimate available 
which originates from the overlap description of domain wall fermions |^5[ (thus it is exact only for 
Ng -^ oo). In this case a transfer matrix can be defined which describes propagation in the fifth 
dimension T = exp(75//^(— M)) where 



H^{-M) = -y^D^-M) 



(40) 



is the ordinary 4(i Hermitian Wilson-Dirac Hamiltonian with a mass term that is negative of 
the conventional one. D" first supports exact zero modes as M approaches a critical value Mc 
defined by a vanishing pion mass. In this case T has a unit eigenvalue and propagation in the 
fifth dimension is unsuppressed. Thus Mc for domain wall fermions corresponds to Mc for Wilson 
fermions, usually given in terms of the hopping parameter, k, 



w. 



M. 



w 



1 



2k 



w 



(41) 



We can simply take Mc = M^ from existing numerical simulations. In Table |^ we give the values 
of k)^ computed in Ref. |36| and the corresponding value of M^ which we use for the rest of this 
work. 

In Figure y we show schematically what the spectrum of the four-dimensional Wilson-Dirac 
operator should look like (with the domain wall convention for the sign of M) for the single flavor 
case. Ordinary Wilson fermion simulations are performed in the region M < Ale- while the domain 



wall fermion simulations have Mi < M < M2. It has been conjectured by Aoki [37| that there is a 
range Mc < M < Mi where there is no gap due to the spontaneous breaking of flavor and parity; 



evidence for the existence of this Aoki phase has been found in lattice QCD simulations |35] and 
analytically using an effective chiral Lagrangian [^9| where the width of this phase. Mi — Mc, was 
found to be 0{a'^). In the conventional picture (see Ref. [39| for a recent summary) the gap reopens 
after Mi and closes again at M2 when four of the doubler modes become massless. 
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M 



Figure 2: Expected spectrum of the Wilson-Dirac operator D" as a function of M (which is 
negative compared to the usual convention for Wilson fermions). There is a mass gap for M < Mc 
and Ml < M < M2, and no gap in the Aoki region Mc < M < Mi. 
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The spectrum of the Hermitian Wilson-Dirac Hamiltonian has been studied in some detail in 



Ref. |4C]. Many level crossings in the Hamiltonian were found uniformly in the region above Mi. 



These zeroes were related to small instanton-like configurations and are presumably related to the 



same lattice artifacts that give rise to so-called exceptional configurations |41]. If the density of 
these zero modes is non-zero in the large volume limit, then the gap is closed and domain wall 
fermions cannot exist: the entire allowed region of M should then be in the Aoki phase. Since 



numerical evidence |E2, 23, 24, ESl for the existence of domain wall fermions is quite strong, we 



infer that the density of these zeroes vanishes in the large volume limit (see Ref. [16| for more 
details), at least for the couplings relevant to present simulations. This is also consistent with 
the above studies of the Aoki phase. However, we also note that at some strong coupling the 
conventional picture is for the whole region to be in the Aoki phase, and at this coupling domain 
wall fermions cease to exist. The existence of the Aoki phase also reveals why the size of the extra 
dimension must increase as the coupling becomes stronger: as the gap gets smaller the size of the 
extra dimension must increase to maintain the same amount of suppression. We refer the reader to 
Ref. [^] for similar plausibility arguments on the behavior of domain fermions at strong and weak 
coupling. 

In closing this section we emphasize that the simple replacement of 

M ^ M = M-MY (42) 

in the quark mass (P6|), though nonperturbative, is an ansatz which may or may not introduce 0{a) 
errors. Furthermore, M cannot really be uniformly shifted, even piecewise, over the whole region. 
Nonperturbative effects such as instanton-like artifacts may be important (though these do seem to 



be more or less uniformly distributed above Mc). However, as previously mentioned, Eqn. (42) is a 
very good fit to our numerical data which we present in Sec. Q. Also the identification in Eqn. ( p|) 
is only exact in the limit Ng -^ cxo. Again, simulations indicate this is a good approximation (see 



Ref. |16|). 



3.3 Quark mass renormalization 

In this Section we concentrate on the renormalization of the quark mass. We follow the method 
outlined in the case of the wavefunction renormalization PQ] . The tree level quark mass was given 



in Sec. 2.2 by finding the smallest eigenvalue of the mass matrix squared, fio^o- ^^ one-loop level. 



the mass matrix is renormalized: 

n = Qo + g^ni. (43) 

i^o is the tree-level mass matrix (|ll]) and J^i is the one-loop correction given by the terms in ( p5|) 

^1 ^ -^(sW.-+amS(2)-). (44) 

where S(j)'± = S(j')p±. 

We rotate the fermion fields to the basis which diagonalizes the one-loop matrix fl 



(p) = Us,s'P+^s'{p)+Vs,s'P-^s'{p) 

^f'Hp) ^ A'{p)P+{V^)s',s + A'{p)P-{U^)s',s- (45) 
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Then the terms which control the renorniahzation of the hght fermion mode are as fohows: 



yiy+C/t 



=l,u=l 



yL(2).+[/t 



1,1 
1,1 
1,1 
1,1 
1,1 



= UW'V^ 
C/l(o),-yt 

y/(i).-Ft 
c/L(2).-Ft 
C//(2),-yt 



1,1 ~ 


mM{2 - M) 


1,1 


0{Nsb^^) - 


1,1 ^ 


u 


1,1 ^ 


h 


T— 1 


M{2 - M) U 


1 1 


M{2-M) h. 



neghgible 



1,1 



The results of Ref. |2C] for Li and Ii combined with our results for L2 and I2 are that 
Li = 



-2 / dx xln( — + - 
io ^azJ 2 

Z2 = -4 / dx Inf— )^ + 4 
^ r dH f 1 V- 



(46) 
(47) 
(48) 
(49) 
(50) 
(51) 

(52) 
(53) 



(2vr)4 ]si2 



sin^ /^(Gij + Gl) + 2 cos /^(6o - &(/))G'ij 



^ sin%, 
. 2Z~4 



^0 - 6(0)Gii - (2 cos^ I - ^ cos^ I) G,. + PGr 



+ 167r^ 



'^ d^/ 



1 



-167r^ 



'^ d'^/ 1 1 



(54) 



_. (27r)4 [2 (1 _ e-"(06o)2 



pj: 



cos 



2 'm 



1 



72 



-a(0 



z,. 



+ ^sin^e-2°(') 



2 ;2 



167r^ 



2 

'^ dH 1 



+ ^<^^'VIilF®<'^-' 



2 (1 - e°(06(/))2 (1 - e"«6(/)) 

1 



(2vr)^ 



l^ 



-. (27r)4 [2 ■ 
,(0)^2 



(55) 
(56) 



where z^ = {l — x){p'^x + {mq )^), /^ = sin/^ and /^ = 2sin(/^/2). We have also used the definitions 



Gr{1) 



2h{l)s].nha{l) 



(6^i-e-"W)-(6o-e°(')) 



l-hl 



and 



Gl{1) 



1 



2&(0sinha(/) 



^i_e-a 



^1 _ e-"(0 ) + (&o - e"(0 ) (e"(0 - 60)2 . 

W)-(6o-e"(')) l-6g e"«-6(0 



*{0) + (60 - e"(0) (e°(0 - 6o)2 e-"W - h{l) 



(57) 



(58) 
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Figure 3: The three finite integrals evaluated at values of M, (see Eqns. (|5^)-(p6D). 



The quantities b{l), bo = b{0), and a{l) are defined in Sec. 2.3 



The finite integrals Ii,l2, and T are 
plotted as functions the five-dimensional mass M in Figure \ 

The expression for the fermion self-energy in the continuum, using dimensional regularization 
in the MS scheme, can be written as 



^MS 



^{p,mf) = ii>T.f'^ + m^Sf^ 



where ruf \s the continuum mass and 



*/ 



j^MS 



vMS 



2 / da; X In 



4/;..ln(|)-2, 



(59) 

(60) 
(61) 



and z^ = (1 — x){p^x + m?f). The lattice mass and continuum mass are matched onto each other 



by combining these calculations so that, in Eqn. (30), 



^ X ^ 



Z\^^ Z. 



MS ' 
2 



(62) 



with the notation that Zj = 1 — (/^Ci^Sj/lGvr^. The final mass renormalization factor between 



16 



u 



4 



3 - 



2 - 



1 - 








0.5 



1.5 



2 



Figure 4: The matching coefficient Cm as a function M. The dashed hues indicate the range of 
M where numerical simulations have been performed. 



domain wall fermions and continuum fermions in the MS scheme is given by 



where 



An 



Un 



1 



luTT 



167r2 



1 1 

4 + 6 



ln(/za) — C„ 



+ h 



(63) 



(64) 



As discussed in detail in Sec. |3.2| , the parameter M is renormalized. Since the perturbative estimate 
of the renormalization ( p8| ) is large and untrustworthy, we use the ansatz ( |4^ ) . The integrals ( p4| ) 
and ( p5|) are unchanged, except the trivial replacement M ^ M in Figs, y and ^ We mark 
the region of M where our simulations have been performed with vertical dashed lines. Between 
those, 0.7 < M < 1.0, Cm varies from 2.88 to 3.10. These values of Cm should be compared to 
Cm = 2.16 for Wilson fermions, 3.22 for Sheikholeslami-Wohlert fermions, and 6.54 for Kogut- 
Susskind fermions. 



4 Monte Carlo results 

In practice, Monte Carlo simulation of quenched QCD using domain wall fermions is very similar 
to standard calculations using Wilson fermions. We use the conjugate gradient algorithm to invert 
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Table 2: Pseudoscalar mass and decay constant for the main numerical data set. 



m 



P = 5.85, 16^ X 32 P = 6.0, 16^ x 32 



/3 = 6.3,243 x60 



0.025 
0.050 
0.075 



0.488(5) 
0.588(4) 



0.104(4) 
0.114(4) 



0.309(6) 0.076(4) 
0.423(5) 0.088(4) 
0.517(5) 0.094(4) 



0.245(5) 0.056(5) 
0.340(4) 0.066(3) 
0.425(4) 0.072(3) 



the five-dimensional fermion matrix. Next, 4(i quark fields are constructed from the 5d fields at 
the two boundaries |18|: 



(x) = P+Mx) + P-^^nAx) 

(x) = '(Aiv,(x)P+ + Mx)P- 



(65) 



These are the simplest interpolating fields for the lightest mode, and composite operators con- 
structed from them satisfy exact continuum-like Ward identities in the limit A'^s — > oo [p^ . 

For this exploratory calculation of the strange quark mass, we compute the pseudoscalar meson 
mass and decay constant on a few dozen configurations at three lattice spacings. Specifically we 
perform Monte Carlo simulations at three values of the gauge coupling: P = 6/g^ = 5.85, 6.0, and 
6.3. The size of the extra dimension for the main part of this work for the three couplings was Ng = 
14, 10, and 10, respectively, and the 5d mass parameter was M = 1.7, 1.7, and 1.5, respectively. 
The value of M ~ 1.7 at 6.0 is the optimal value which suppresses propagation of the light mode 
in the extra dimension, as found in a previous study p3| ]. The M at the other two /3's were ad hoc 
choices based on the P = 6.0 value and the fact that the optimal M should decrease to 1 in the 
weak coupling limit. At all three gauge couplings the mesonic two-point functions were computed 
with m = 0.075 and 0.050, and at /3 = 6.0 and 6.3 another mass, m = 0.025, was also included. 
The lattice volumes at /3 = 6.0 and 6.3 are roughly ~(1.6 fm)^ while the P = 5.85 volume is ~(2.0 
fm)^. Our raw lattice simulation results are given in Table |2[ 

In addition to the gauge coupling and bare quark mass, domain wall fermion simulations depend 
on the five-dimensional mass M and the size of the extra dimension Ng. The study of how data 
are affected by changing these parameters is important. Given the results of Ref. ||2^, we believe 
that we have taken A^^ large enough so that corrections due to the mixing of the two light modes 
in the center of the extra dimension are smaller than the rest of our uncertainties which we will 
show later to be ~ 10% — 20%. The effect of a finite Ng is to give the pion an additional intrinsic 
mass mi: 

{aM^f = AM {2-M) {m + mi). (66) 

In the free theory mi = |1 — M|^ [|l9[, and in the interacting theory m/ is expected to decay 
exponentially m/ ~ exp(— aA^^). 

Indeed we find a nonzero intercept when we extrapolate the pion mass squared to m = for 
P = 5.85 and 6.0, (see Fig. ^, but find that the P = 6.3 mass squared does extrapolate to zero. 
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Figure 5: Lattice pseudoscalar mass vs. m for the three f3 values. Asterisks indicate Unear 
extrapolations to m = with jackknifed errors and are slightly displaced horizontally for clarity. 



We have repeated simulations at /3 = 6.0 using Ng = 14 and found the pion mass decreases by 
a few percent (compare Tables § and y) signaling a statistically significant mj. However, errors 
in determining the lattice spacing from /^r are much larger than the difference in aM^^ and so the 



N. 



10 results are sufficient for this work. 



A thorough study of the Ng dependence at a stronger coupling, (3 = 5.7, was presented in 
Ref. 1 24, 25]. They find two important results. First, physical results were unchanged between 
Ns = 32 and 48. This is consistent with arguments above concerning the behavior of domain wall 
fermions at stronger coupling. They used M = 1.65 in their study. It is possible that a larger value 
would decrease the value of Ng required to reach the asymptotic region. Second, even in this limit, 
the pion mass does not vanish as m — > 0, which is then presumably a (4d) finite volume effect. 

We were able to test for finite volume errors at /3 = 6.0 by computing the pion mass with 
Ns = 10 on a 24^ x 40 lattice. The results, aM^ = 0.318(5) and 0.427(4), for m = 0.025 and 
0.050 respectively, are not statistically distinguishable from the results on the 16'^ x 32 lattice (see 
Tab. 1). 

Since the quark mass, even at tree level, explicitly depends on M, we performed several more 
runs at /3 = 6.0 with m = 0.025,0.050 and 0.075, varying M between 1.5 and 2.1 on roughly 20 
configurations with Ns = 14.Q Table ^ displays the values of aM.,^ obtained from these simulations, 
and Fig. ^ shows (aM-,^)^ as a function of m. The linear least squares fits extrapolate to (aM-,^)'^ = 

^ In the interest of frugality, we did not extend the m = 0.075, M = 1.7 data from Ns = 10 to Ns = 14. 
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Table 3: Pseudoscalar mass (aM,r) for the (3 = 6.0, Ng = 14 numerical data set. 



m 



M = 1.5 M = 1.7 M = 1.9 M = 2.1 



0.025 
0.050 
0.075 



0.290(9) 
0.394(7) 
0.487(7) 



0.293(7) 
0.411(5) 



0.290(7) 
0.412(7) 
0.510(7) 



0.281(11) 
0.396(8) 
0.490(7) 



at ?7i = within the (uncorrelated) errors for each M separately. Furthermore, in Fig. |^ we plot 
{aM-j^y' as a function of M for the three values of m. We can test the ansatz (42) simply by fitting 
the data to 

{aM^f = AmM{2-M), (67) 

for each m. The dashed lines in Fig. |^ are fits to Eqn. ( |67| ) and have good x^'s; each x^ P^r 
degree of freedom is less than one, and A is the same within (uncorrelated) errors for each m: A = 
3.48(8), 3.42(6), and 3.52(4) for m = 0.025, 0.05, and 0.075, respectively. Therefore, the ansatz 
( [4^ is well justified in this work. At present the data do not permit a more general fit. Given 
these observations we find the definition of the lattice quark mass 



m 



LAT 



m M (2 - M) 



(68) 



to be very reasonable and suitable for this work. 

5 Quark mass and coupling constant 

We believe we have a good understanding of the dependence of hadron spectrum on M, so now we 
concentrate on the data at one value of M per lattice spacing (as listed in Table |^. We compute 
the pion decay constant afj^ for three masses m for each /? and extrapolate linearly to the chiral 
limit, 171 = 0, (see Fig. ^). We determine the lattice spacing by setting this extrapolated value to 
the physical pion decay constant, /,r = 130.7 MeV. Here we must emphasize that in this exploratory 
calculation, the determination of the inverse lattice spacing has a large systematic uncertainty due 
to having a only a few data points to extrapolate af-jt to m = 0. In the continuum, xPT gives a 



one-loop correction to the decay constant which goes as [42| 



f. = f 



1 + 



Mj 



(4vr/)- 



■In 



A2 



(69) 



where / is the tree-level decay constant and A is the cutoff. With only three quark masses we 
cannot resolve the logarithmic behavior of /^r, so we extrapolate linearly to ?7i = 0, (see Fig. |8|). 
We should remark that the signal for the decay constant is rather noisy: a/,r varies by as much 
as 10% depending on the range of Euclidean time over which one computes correlation functions. 
Since the determination of the lattice spacing for this work is through the extrapolation of af-,t to 
m = 0, it comes with a large uncertainty. 
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Figure 6: Pion meson mass squared as a function of ra for /3 = 6.0. Lines are least squares fits 
which extrapolate to (aMyr)^ = at m = for all four values of M within the errors shown by 
asterisks. 



Table 4: Summary of simulation parameters and results. 





13 = 5.85 


/3 = 6.0 


13 = 6.3 


# configs. 


18 


30 


11 


volume 


16^ X 32 


16^ X 32 


24^ X 60 


iV, 


14 


10 


10 


M 


1.7 


1.7 


1.5 


a-^U) (GeV) 


1.57(15) 


1.89(14) 


2.72(34) 


m,(MeV) 


73(7) 


75(6) 


76(10) 


Mf 


0.908 


0.819 


0.708 


^m 


2.94 


3.01 


2.94 


mLAT(MeV) 


70(7) 


74(6) 


73(10) 


(Tr t/piaq/3) 


0.5751 


0.5937 


0.6224 


aviSAl/a) 


0.157 


0.146 


0.131 


Z^ili = 2GeV) 


1.30 


1.34 


1.34 
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Figure 7: Pion meson mass squared as a function of M at /? = 6.0. Ng = 14 except for the M = 1.7, 
m = 0.075 point which has Ng = 10. The dotted hues are fits to (aMT,)'^ = AmM{2 - M). 



Next, we use chiral perturbation theory as a guide to interpolate the pion mass squared hnearly 
in the quark mass: M^ ~ m to the value of m which gives the physical kaon mass, 495 MeV. That 
value must then be doubled to take into account the fact that we use degenerate quarks in our 
simulations, while the physical kaon has one strange quark and one lighter non-strange quarkj^ 
We denote the parameter m corresponding to the physical kaon by rUg and display our results in 
Table ^. Using Eqn. (pq ) we obtain m^ , the domain wall strange quark mass. 

Combining Eqns. ( |30[) and (|6^) the following expression relates a quark mass computed on a 
lattice with spacing a to the quark mass defined in the modified minimal subtraction scheme of 
dimensional regularization at momentum scale //: 



m 



MS 



(/^) 



m 



LAT 



1 ( ln(/ia) — C„ 

TT 



(70) 



where we have substituted 5^/47r 



Q. 



and Cp = 4/3. The last quantity we need is the coupling 



constant a^. In a one-loop calculation such as this, there is ambiguity in the definition of the 
coupling constant. The matching equation ( [70[ ) is derived by equating the poles of the one-loop 
quark propagators computed in the continuum and on the lattice. Each procedure uses a differently 
defined coupling constant; however, the difference between the two coupling constants is a higher 

^To be accurate one should use m = irna + mi)/2, setting mi with the physical Mtt. However, in this exploratory 
work, we have neglected this ~ 4% effect. 



22 



0.15 



0.1 



a 



0.05 








D 


5.85 


O 


6.0 


o 


6.3 







0.1 



0.05 

m 
Figure 8: Lattice pseudoscalar decay constant vs. m for the three (3 values. Asterisks indicate 
Unear extrapolations to m = with jackknifed errors. 



order correction in perturbation theory: 



a 



CONT 



0{a 



(71) 



It has been known for some time that the bare coupling constant is not a good expansion 
parameter for lattice perturbation theory |^^ . Therefore we use a physical definition related to the 
heavy quark potential. Specifically, we define the coupling constant using the two-loop perturbative 
expression for the plaquette (the 1x1 Wilson loop): 



In ( -Tr a 



plaq 



47r 



ay(3.41/a) 1 - 1.19av + 0{a 



(72) 



The scale 3.41/o has been computed by estimating and minimizing the effect of higher order 
terms p5| |. One can run ay to any other scale using the universal two-loop beta function. Al- 
ternatively, one can compute the continuum MS coupling constant a^jg from ay (at a scale q) 
perturbatively E^]: 



«Ms(9e ^^^) = OLv{q) 



1 



2av 



vr 



+ 0{a 



Although the three-loop beta function has been computed for arro 



(73) 
we choose to use the 



two-loop beta function for consistency; the difference comes in below the 1% level. 

The remaining problem is to decide on the scales and corresponding couplings to insert in the 
matching equation (^). Refs. [^, 4£] advocate the reorganizing of lattice perturbation theory 
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Table 5: Strange quark masses (in MeV) in the MS scheme at increasing values of the matching 
scale fi. We give the mass at both the matching scale fi and at 2 GeV. The first pair of masses are 
without tadpole improvement and the second pair include tadpole improvement. 



f^ 



aj-r^ifi) m{n) m(2 GeV) m^^fi) m^\2 GeV) 



,TI/ 



IMS 







/3 = 


= 5.85 






0.5/a 


0.275 


114(11) 


96(9) 


108(11) 


90(9) 


l/a 


0.199 


96(9) 


92(9) 


93(9) 


90(9) 


2.0 GeV 


0.181 


91(9) 


91(9) 


90(9) 


90(9) 


2/a 


0.156 


85(8) 


90(9) 


85(8) 


90(9) 


vr/a 


0.137 


81(8) 


90(9) 


81(8) 


90(9) 







P 


= 6.0 






0.5/a 


0.243 


116(9) 


102(8) 


111(8) 


97(7) 


l/a 


0.182 


99(8) 


99(8) 


97(7) 


96(7) 


2.0 GeV 


0.178 


98(8) 


98(8) 


96(7) 


96(7) 


2/a 


0.146 


89(7) 


97(7) 


89(7) 


97(7) 


n/a 


0.129 


85(6) 


96(7) 


85(7) 


97(7) 







/? 


= 6.3 






0.5/a 


0.202 


107(14) 


101(13) 


103(14) 


97(13) 


2.0 GeV 


0.175 


99(13) 


99(13) 


97(13) 


97(13) 


l/a 


0.159 


94(12) 


98(13) 


93(12) 


96(13) 


2/a 


0.131 


86(11) 


97(13) 


86(11) 


97(13) 


n/a 


0.118 


82(11) 


97(13) 


83(11) 


97(13) 



as we described above, so that the resulting expression may be "horizontally matched" to the 
continuum perturbative expansion. In converting our lattice quark mass into a continuum quark 
mass we follow the procedure described in Ref. ^. Then the continuum matching scale /x should 
be set to the "best" lattice scale which minimizes the higher order corrections to the fermion self- 
energy. Unfortunately, it is harder to estimate ^ for logarithmically divergent graphs than it was 
in the case of the plaquette. Therefore we resort to trying a spread of /x values. Evidence from 
previous work indicates that in the range ^ = 0.5/a to vr/a the higher order, ultraviolet-dominated, 
effects are minimized |45, ^. Finally, the quark mass is run to 2 GeV using the two-loop running 
equation |Q. We also test the systematic error by repeating the procedure for different values of 
fi. In Table ^ we give the strange quark masses at the matching scales, m{fi), and the mass run to 
yu = 2 GeV; we also give our values of a^jg at different scales. Our statistical errors are 3-4 times 
larger than the variation in m{2 GeV) computed with different values of //, so we cannot argue that 
one particular scale minimizes higher order effects in Eqn. (|70|). 

One might try to improve the perturbation expansion by nonperturbatively estimating ultra- 
violet effects which spoil the convergence. For example, the tadpole-improvement prescription 
advocates perturbatively expanding the quantity mr /uq, where uq is designed to be sensitive to 
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short-distance fluctuations [^]. Therefore, the tadpole-improved quark mass, in the MS scheme 
is related to m^^'^ /uq perturbatively by 



m^if^) = {m'^^^/uo) 



1 ^"« I 1 

1 m fia 

vr 



Cm — Ct 



(74) 



where uq is defined to be the fourth root of the plaquette 



/ 1 \ 1/4 

UO = ( oTr ^plaq) , (75) 

and it is computed from Monte Carlo simulation. The matching coefficient is modified at the one 
loop level by the perturbative expansion of uq 

1 2 

uo = 1 - —-g^ = I CtUs, (76) 

II vr 

which defines Ct- In Table H we give the tadpole-improved strange quark mass for the various 
matching scales, rn^^{^), as well as the mass run to 2 GeV, 171^^(2 GeV). At the lower scales /U < 2 
GeV, tadpole-improvement lowers the mass by 4-6 MeV indicating that the unimproved pertur- 
bative result has significant higher order corrections. On the other hand, there is no significant 
difference between the standard and improved expansions when the matching is done at ^u > 2 
GeV. Therefore, we will choose for our final result the improved mass with the matching done at 
/i = 2/a, m^(2/a) in Table |5[ and assign a 2 MeV systematic error due to the arbitrary choice of 
scale. We perturbatively run our final result to 2 GeV for comparison with other results. 

In Figure |9| we plot our results for the strange quark mass (with statistical error bars only) in 
the MS scheme at the scale 2 GeV along with those obtained using other fermion discretizations. In 
choosing the data to which we compare the domain wall results, we used those which were obtained 
with the same method. The masses were computed by fixing the kaon to be its physical mass and 
the matchings to the continuum were computed perturbatively. Unfortunately, we can only set the 
lattice spacing using f^^ while the other data in Fig. ^ use the p meson mass. In typical Wilson 
fermion simulations [51|, as well as preliminary domain wall fermion simulations |^^, there is a 
?a 15% systematic uncertainty due to choosing 7,^ vs. Mp to set the scale. 

In this first work, the statistical errors are rather large, about 15%, yet it is encouraging that 



we see no significant signs of scaling violations. As argued in Sec. p.3| , discretization errors should 
be 0{a^) rather than 0(a) since wrong-chirality operators are suppressed; however, a more precise 
study is needed to draw a firm conclusion. Since no scale dependence can be detected within the 
statistical errors in this exploratory study, we take a weighted average of the strange quark mass 
determined at the three lattice spacings to give 95 MeV with a purely statistical uncertainty of 
5 MeV. A linear extrapolation to the continuum limit raises the central value by 15 MeV, but of 
course has a large uncertainty. Given that we do not expect to have 0(a) scaling violations, we 
take 95 MeV as our strange quark mass but add the 15 MeV in quadrature with the rest of our 
systematic errors. 

The largest systematic error arises from the ambiguity in which physical quantity is used to fix 
the lattice spacing. As mentioned above, preliminary results and experience with Wilson fermions 
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Figure 9: Strange quark mass in MS scheme at 2 GeV. Our results are displayed as squares, Wilson 
fermions as diamonds [|^], Slieiklioleslami-Wohlert fermions as circles [^, and Kogut-Susskind 
results as crosses |43, ^]. 
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Table 6: Sources of systematic uncertainty in computing the strange quark mass. 



Source 


Estimated size 


Using /tt vs. Mp to set a 


20% 


Matching scale, pL 


2% 


Using (j) vs. K to set rng 


10% 


Continuum extrapolation 


16% 



lead us to believe there could be at least a 15% uncertainty using domain wall fermions. Considering 
that our calculation of f-^ was done with only a few quark masses, we conservatively estimate a 20% 
systematic uncertainty in a~^; the other systematic errors are small in comparison. For example, 
the uncertainty in matching scale induces an error of 2 MeV. We have not yet computed the strange 
quark mass by fixing the (j) meson to its physical mass, but in Wilson fermion simulations there is 
an 0(10%) difference from the strange quark mass using the kaon. For now we include this error 
as a systematic uncertainty in our calculation, but it must be explicitly checked with domain wall 
fermions. Of course quenching also induces a systematic error; however, we are unable to address 
its effects without a full dynamical-fermion simulation. 

Adding the statistical and systematic errors in quadrature, our final result for the strange quark 
mass within the quenched approximation is 

mj^(2 GeV) = 95(26) MeV. (77) 

6 Conclusions 

This work has focused on the steps needed to compute the light quark masses with the domain wall 
fermion discretization. We have extended the calculation of the domain wall fermion self-energy 
to the massive case. We find that the perturbative mass renormalization factor which matches the 
domain wall lattice regularization to the MS regularization scheme is as well-behaved as that for 
Wilson fermions. 

In conjunction with the perturbative calculation, we have performed numerical simulations of 
quenched lattice QCD using domain wall fermions. At /3 = 6.0 the pion mass squared vanishes 
linearly in m as m ^ 0, and the ansatz M ^ M — M^ is a good fit to the data. Finally we 
compute the strange quark mass my^{2 GeV) at three lattice spacings. Within our errors, the 
results are scale independent so we take a weighted average giving r7T,j^(2 GeV) = 95(26) MeV, 
where systematic uncertainties (except for quenching effects) have been added in quadrature with 
the statistical error. 

We intend to perform a larger scale calculation which will include the vector meson spectrum. 
This will allow us to calculate the "average" light quark mass {niu + m^jl together with the 
strange quark mass. In addition we will be able to estimate the systematic error due to setting the 
lattice scale from /^ vs. Mp. A higher statistics study will be able to sensitively test for scaling 
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violations, and possibly give a value for the light quark masses which is comparable in precision to 
other lattice calculations. 

Future prospects for this formulation are promising. Domain wall fermions have an advantage 
over Wilson fermions, improved or not, in that they have chiral symmetry which is broken only 
by the explicit mass m coupling the boundaries - the mixing of the two modes in between the 
boundaries can be made negligible compared to m. Consequently there is no mixing between 
operators of different chirality and the quark mass is protected from additive renormalization. 
A lattice discretization which respects the axial symmetries of continuum QCD has an excellent 
chance to improve calculations of matrix elements involving light hadrons. 
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A Diagonalization 

A.l Tree level 

In this Appendix we discuss the spectrum of the domain wall Dirac operator at zero momentum 



at tree level. Although this work appears in Ref. |17|, a treatment here is useful to establish our 
notation, which is similar to Ref. |20|. 



Given the mass matrix i7o, defined in Eqn. (|ll|), we wish to solve the eigenvalue equation 

(n^nl)^/:} = (A«)V«, (78) 

where the index i € [1, A'^s] labels the eigenstates. The general equation to be satisfied is 

(i + 62_a2)0, -6o(0,+i+0,_i) = 0. (79) 

If < 6o < 1 then the general solutions are of the form 

cj), = Ae"»^ + Se-°o' (80) 

where ao is defined through 

If — 1 < 6o < then the general solutions are of the form 

(Ps = (-l)"(^e°o" + 5e-"''^). (82) 
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For (Aq )^ < (1 — 60)^ or (Aq )^ > (1 + bo)'^, a becomes imaginary and 



W^2 



COS laQ 



l + bl-j X'o') 
2bo 



(83) 



Let us concentrate on the first case of exponential damping. We must apply the boundary 
conditions 



(bl + a^m^ - (A5^V)</'i - bo 4>2 + am bo (pNs 
[l + bl- {X^fjct^Ns - bo 4>Ns-i + am bo (pi 




0. 



(84) 
(85) 



In order to make these conditions consistent with the general equation ([79|), 08=0 and 4>s=Ns+i 
must satisfy 



^0 00 + 01 
-bo 4>Ns+l 



(am,) 01 — am, bo ipNs 
am, bo 4>i- 



(86) 

(87) 



As has been found previously |T^, ^, IS], the eigenvalue is 

(A«)2 = {am)\l-blf + 0[{amf) + o[b^') 

= {amfM^{2-Mf + o({amf'^ + o((l-M) 

The corresponding normalized eigenvector is 



Ns 



^'P 



M{2-M) e-"«("-i). 



(89) 



The eigenvectors corresponding to the heavier modes can be decomposed into a basis of sine func- 
tions 



di) 



Note that if 60 < then 




2 /7r(i-l)r 



Ns + l-s 



; for i ^ 1 



0W = (-i)«yrT^e-"»(^-i) 



(90) 



(91) 



while the other eigenvectors and all the eigenvalues are unchanged. Also, the eigenvalues and 
eigenvectors are the same for i^Qr^o with the substitution s ^ Ng + 1 — s. Let us define unitary 
matrices which diagonalize Oof^Q as 

(92) 

(93) 






/,(«) 



and which diagonalize rijilo as 



s,s' 
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A. 2 One— loop level 



As in Appendix |A.l| , we want to derive an effective action for the lightest eigenstate. In general 
this involves computing one-loop corrections to the matrices U and V which diagonalize the square 



of the mass matrix, $7. However, it has already been shown [^, 21] that the light mode is stable 
under radiative corrections. Following Ref. [Eol, let us write the mass matrix as 



n 






(94) 



where the T.^^^ are given in Eqns. (32) - ( ^4|) and S^-')'^ = T,^^')P±. Then, to diagonalize OO"!" to 
0{g^), one must compute the corrections to U and V: 



[/(O) ^ [/^(l+^2^(l))[/(0) 
y(0) ^ y=(l+^2y(l))y(0)^ 



If we write the one-loop eigenvalue equations as 



(uQ,n^u^\ 



J s.s' 



{\i'^f+9\x['^r 



and 



ixi'^f+M'^f 



V J s,s' 

then it can be shown that ||2^ 

"(c/(°)f]iyWt),,,,A(^') + (f/Wf^iyWt),, ,a(^) 



(95) 

(96) 

(97) 









(aW)2 - (A^V 



for s y^ s', and 

and, most importantly. 






0, 



(^W)2 ^ 2([/Wj7iyWt),,,A(^\ 



(98) 



(99) 



Therefore, in order to compute the quark mass to one loop, we need {U^^'QiV^^'^)i^i. Aoki and 

Taniguchi pO] have shown that U^^'Tjq^^V^^'' is negligibly small, and in Appendix O we compute 

[7(0)s(2).-y(0)t andy(0)s(2).+C/(0)t' . 

1,1 



1,1 



B Perturbative calculation — details 



B.l Feynman rules 



The Feynman rules are similar to those for Wilson fermions plus plaquette-action gluons [53| 
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In the gauge sector, we use the usual four dimensional plaquette action: 

Sg = /3E^I^eTr(l-C/pi,q). 

plaq 



(100) 



In the five dimensional picture there are simply Ng copies of the gauge field, and in the flavor 
interpretation the gluons are simply flavorless. 

The propagator of a gluon (in Feynman gauge) with momentum q is given by 



D^^uiq) 



JflU 



4EpSin^(a(?p/2) 



-'fJ.U 



(101) 



The one-gluon - fermion vertex with incoming (outgoing) momentum -pi (^2) is gtnVfiiQ), where 

Vfiiq) = (ijf,cosaqfj,-sinaq/j, (102) 



Q = (Pi +P2)/2, and t„ is one of the eight generators of SU{3). The two-gluon - fermion vertex 
with gluon indices n,fi and n, fi' and fermion momenta as above is (5^/2){t„,f„'}(^^^p', with 



Wf,iq) 



(^ijfj, sin aqfj, + cos aq^j. 



(103) 



The massive fermion propagator for the boundary wall variant used in this work was derived 
in Re f. p^ and also Ref. [^0| whose notation we adopt. The tree-level G and G , as defined in 
Eqn. (|20|) , are given by 

Gf ,, = gL' + A++e"(^+^') + 2yl+_ cosh a{s - s') + A__e-"(^+^') 



and 

G^^^, = G°3, + S++e"('+'') + 2S+_cosha(s-s') + ^— e^ 
The inhomogeneous part of the solution is 

Gg^gi = Acosh.(a{Ns — \s — s 

and the constants defined implicitly above are 

A 



-a{s+s') 



A)); 



1 



A+. 
B+- 



46sinha sinh(aA^s) 

26 sinh q(1 + 2m cosh(aA''s) + m'^)A/F 
(1-e 

(g-2a7V, 

26 sinh a(l + 2m cosh(QA''s) + m^)A/F 
(1 _ e2"^^)e"(e" - 6)(1 - m^)A/F 



^"^^){l-be''){l-m^)A/F 

l)e-"(e-"-6)(l-m2)^/F 



where 
F 



^aNs 



1 -be"' + 711"^ {be-'' -1) 



4:bm sinh a + e 



-aNs 



be~ 



l + m^il-be° 



(104) 
(105) 

(106) 

(107) 

(108) 
(109) 
(110) 
(111) 
(112) 
(113) 

(114) 
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B.2 Tadpole diagram 

The tadpole diagram (Fig. |l|a) is simple to compute. Since it has no fermion propagator in the 
loop, it is equivalent to the case for Wilson fermions: 



-(¥+-)t, (115) 



where T is the finite integral 

T = 167r2 / -^ ,, , , = 24.4. (116) 

y_. (27r)4 4E;.sin2(/^/2) ^ > 

B.3 Taylor expansion 

The calculation of the half-circle graph (Fig. |l|b) with domain wall fermions parallels the same 
calculation with Wilson fermions [34|. The contribution of the half-circle graph to the fermion 



self-energy is 

S^7(P,-) = ^/_'^^(^E-M(^^)C'(U-)I)..(ap-0^'.(^^). (117) 

-E{l,ap,am). (118) 



li,u 

167r2 r dH 



a J-T, (27r) 

The second line above defines the integrand E. 

A simple Taylor expansion of Eqn. ( |117| ) about zero lattice spacing would not be valid due to 
the logarithmic divergence of the integral. That is, the coefficients of the power series in a would 
have a In a dependence, which must first be separated before expanding the Taylor series. We 
subtract and then add back a similar half-circle graph built from continuum-like Feynman rules 
designed to have the same infrared behavior as the present rules. 

.IR . ^ 16^2^- dH ^, ^ ,,^ (^7^) 5i^,(/,am) (i7^ 



-E^^(l,ap,am). (120) 



a y-^(27r)4' 

The step function Q{t{^ ~ l'^) niakes the integral spherically symmetric and therefore easier to 
evaluate. We will specify the continuum-like fermion propagator S in the next Section. 
Now we can write the half-circle graph in terms of a Taylor expansion of (S — S ) 

levr- 



a Ji 
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1 

a Ji 
levr- 



E{1, ap, am) - 6(71^ - r)E''^{l, ap, am) 



G(7r^ - l^) E'^^il, ap, am) - E^"(/, 0, 0) + - / E{1, 0, 0) 
a Ji L i a Ji 



da 



E{1, ap, am) - e(7r^ - l^)E''^{l, ap, am) 



a=0 



+ 0(a). 



:i2ii 



We use a shorthand notation for the loop integral, namely /^ = J dl/ {27r). Eqn. ( 121 ) is the 
"master" equation for computing the half-circle graph: in the remainder of this Appendix, we 
compute the various terms appearing within it. 



B.4 Infrared terms 

In this Section, we compute E^^ and dE^^/da. The IR (a — > 0) limit of the fermion propagator is 

1 ^. .., . . .„ „ ^^22) 



5^"(/,am) 



/2 + a2(m^°))2 
where 6+ = ds,ids',N,, S- = ds,Njs',i, A± = dgz^zi^s' - boSs,s', and 



Y,i-if-am6± + A±)C±P±, 



C+ = {l-ht)b, 



2\i2Ns-s-s' 



C- 






(123) 
(124) 



,(0) 



Note also the presence of the tree-level quark mass mq = m,M{2 — M) which was defined above 

in Eqn. {^. 

Let us first compute the integral of dE^^ /da. Since the IR vertices are just i^^, there will be 

just two terms, one with dS^^{l, am) /da and the other with d{ap — l)~'^/da. The result for the first 

term is 

(z7^)5^^(/,am)(-2p) • {ap - /)(i7^ 



©E 



{ap - ly 



a=0 



©E 



-ii>clP± 



l^ 



(125) 



where C^ = C\a=Q, and Q = Q{'k'^ ~ l'^)- Note the A± term in the fermion propagator is anti- 
symmetric in / and so vanishes upon integration above. The second term comes from taking the 
derivative of S^^: 



d_ 
da 



S''^{l,am) 



a=0 



5^^(/,0) 



^(S'''{l,am" ' 
da V 



a=0 



5^^(/,0) 



Y: ^ Hf+ A±) CiP±5± {-i/Ci + AiCO )P± 



(126) 



After some algebra, we have 



(nM)^5'^(^7M) 



a=0 



-4(16.2)^(0) IqI_3 



),N.~s+s'-lp^ 



,Ns+s-s'-l 



P 



(127) 
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The integration of the E {l,ap,am) — i? (/,0, 0) term of Eqn. ( |121| ) is quite similar to that 
for Wilson fermions (see e.g. ||3J]). We delay writing down the answer, since it simplifies greatly 
upon diagonalization in flavor space. 

B.5 Finite terms 

Let's first look at the numerator of E{1, ap, am). For the time being we suppress the indices s, s'. 

= Y. K ''°' ^(^ + ''P^'' ~ '^"^ ^(^ + ''P^^'] [("'^+ W-)G'^P+ + {-i}+ W+)G''PJ 

X |^z7^cos-(/ + ap)^ -sin-(^ + ap)^ (128) 

Multiplying the factors of TV gives 



M 



- E(-^7m?7m + W-) G^P. (cos2i(Z + ap)^ 

- E(-^7m?7^ + W+) G'^P^ (cos2i(/ + ap)^ 
+ ^-ij + W) G^'P^ ( ^ sin2 U1 + ap)^ 
+ (^.ff + w+) G'^P^ (^sin2i(/ + ap)^ 

- E ( - ^^ + ^~)^7/. G^^- {\ Ml + ap) 



- E(-^'^+ ^' 



^7/. G P+ (^-sin(/ + ap)^.y, 



(129) 



where, as usual, 1/ = Yliulylv To compute /;£'(/, 0,0) we divide Eqn. ( |129 ) by r and integrate. 
Since the integration region is symmetric in /^, the terms odd in l^ in ( |129| ) vanish upon integration. 
The result is given by 

i?(/,0,0) = ^i|^sin2|[(M^o-G^),,,,P+ + (t^o'G'^W^- 

- Y.^os'kVWoG'')s,s'P+ + {W,G'')s^,p]-^iGl, + Gl,)\ (130) 
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where W^ = W^\am=o- 
Next we compute 



-—E{l,ap,am) 
da 



a=0 



where 



i]^S^ = Y.'^v^{{l + ap)/2)S^{l,am) — [v^{{l + ap)/2)yi-ap) 



:^ + me"^, 



Ir 



da V 

dS^ 



v^{{l + ap)/2) -^v^{{l + ap)/2) {I - ap) 



and mS'^ 
Note that v^S {dv^/da) is the same as Eqn. (129) with the replacements 



cos^ -(/ + ap)^ 

siri^ -{l + ap)^ 

- sin(/ + ap)^ 



— ^sin(/ + ap)^ 
-^ sin(/ + ap)^ 
-^cos{l + ap)^. 



After some algebra, the result is Sl ^, = 8^ \, P+ + 8^ \, P_ , where 






JX'^^i'^ls' + Gl,) + ^ cos/^((VFo-G^),,,, + {W+G' 



+ 



2/4 



12 



E--'t)^?«'-E^.(2-«'|-E-o«'|)g: 



+ -[{W^G^)s,s' + {W^G'^)s,s 



and 



/<•; = -l[-i^Y{Gl, + Gl,)+Y^co,l,{{W^G^)s,, + {W^G' 



2h 



Ifj, 



'""(E-- Y)G'.'-E^M(2-s^f-E-«'|)Gf 



+ UiW,-G%,,, + {W,+G^), 



:i3i) 



(132) 



(133) 



(134) 



(135) 



For the second term in Eqn. ( |132D , the calculation is simplified through the use of the identity 



-—5, ,'(/, am) 
da '^ ^ 



-Ss,s"{l,0) 



a=0 



-—Dg" s"'(/, am) 
da 



a=0 



Ss"',s'{hO) 



Ss,s"il,0) 1^(63", NsSs"',lP+ + S, 



s",l<Js"'N 



p_ 



Ss"',s'{l,0). (136) 
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The term generated by taking this derivative is proportional, therefore, to the mass m. We express 
the coefficient by 



fj. ' 



v,[^). 



Algebraic manipulation similar to above gives S^ ^i = f^'^, P+ + S^'^iP- with 
Ifi^ = /^{PE-«'|G«^.Gt, -/--Esin^|Gf,iG^,,, 



and 



P 
2 



G^,nAW,-G%^s' + iW,-G%,NA + G^AW+G' 



)Ns,s' 



+ {W+G'^UiGl 



Ii^-^ = /^{pE-«'|Gf,,G«^,, -PE--'|G^^.Gt, 



+ Y.'m'fiW+G%,^{W+G' 



T2 






(137) 



(138) 



Gf,NSWoG%,s' + iW,-G%,N.Gf^,, + G^^AW^G'^U^,, 



+ {W,+G%,,G^^^,, 



(139) 



B.6 Diagonalization of the self— energy 

Returning to the calculation of E^^{l,ap,am) — £'^^(/,0,0), we diagonalize the IR fermion prop- 
agator. We refer the reader to Appendix lAl for the explanation of the matrices which diagonalize 
the mass matrix. In the a — > limit, the propagator for light mode is 



^IR 



S\\{l,am) 



1 



E(-«f+am), 



since the pieces of the full bd propagator (|122|) transform as follows in the diagonal basis: 



UG+U^ = VC-V^ = l + 0{Nsb^') 
VA-G-U^ = UA+C+V^ = 0{N^bl^') 
V6-G-U^ = U6+G+V^ = bo{l-bl)am. 



(140) 



(141) 
(142) 
(143) 
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Then the calculation of the -E'ir(/, ap, am) — £'ir(Z, 0, 0) term proceeds identically to that in the case 
of ordinary four-dimensional Wilson fermions [p3^, pj] . The result is 



levr^ 



a Ji 



ein'-l') Ei^{l,ap,am) - Sir(/,0,0) 



dx 







X <^ li-dx 



azJ 2 



+ 4m^°) 



ln(^)'-l 



i^Li + m(°)L2 



(144) 



(0)^2^ 



where z^ = {1 — x){p^x + (m. 

Next we compute the derivative term of Eqn. ( |121| ) for the light mode. Recalling the notation 
from Eqn. (|33|), 



16vr2 / A 



I da 



a=0 



E{1, ap, am) - G(7r2 - l'^)E^^{l, ap, am 

Combining £^ (Eqns. (|134| ) and ( |135| )) and its IR counterpart ( |125| ) gives 

dH i 1 



ipli + m^"'l2. 



.1°)/. 



/i 



167r" 



^ (27r)4 I 8/2 



E 



sin^ Z^(Gij + Gl) + 2 cos /^(6o - 6(0)^/? 



^ sin%, 
^ 2/4 



(6o - bil))Gn - (2 cos2 | - ^ cos^ |) G^ + PG^^ 



167r" 



'^ r/4/ 1 



where 



^.^(0 - E(f^^°^)MG'5,(0(^^°» 



;s',i 



T.iy'%,sGi,mv^'^^ 



26(/)sinha(/) 



(^, 



-1 _p-"{Oi 



{bo - a 



cc{lY 



l-hl 



{bo^ - e-°(0) + {bo - e"W) (e"(') - 6o)^ 



and 



Gl{1) ^ ^(C/(o))i,,Gf,,,(0(^^°" 



)s',l 



T.(y^'\sGf:Ai){v^'^')s',i 



1 



26(/)sinha(/) 



(feo"' 



'"(0 



)-(6o-e"«) 



1-6^ 



= 0(0 



6(/) 



(6(^^ - e-"(0) + (6o - e"(0) (e°(') - 60)^ e-"« - b{l) 



Likewise, combining S'^ (Eqns. (|138| ) and ( |139D ) and its IR counterpart ( |127D gives 



1 - 62 y_, (2^)4 2/2 



(145) 



(146) 



(147) 



(148) 



''"' ''^ ^'^ :l-{-PT.^os'T^^^sGf,NAuK'^ + U,,GlGl^,vK',i 
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P^sin^f Vi,sG^^,Gf,^^,,U^,,,N^ + U,,sGf^^Gl,,V^ 



s'Ns 





t' 


+ 






2 




P 


+ 






2 




P 


+ 






2 




F 


+ 


— 




2 



y^Gf^N. (W~G^)^^UK',1 + U,,sGf^^^ {W-G^)^^y^ 



s'l 



s',1 



V,,sGl{w+G^)^^^^uK,^ + U,,,Gl{w^G^)^^^^yK,, 
Vi,s{w+G'^)^G§^^,,UK,,, + Ui,s{w+G'^)^G§^^,,V^s',i 

^-^ 2 V /s,l\ ^ Ns,s' 



E 



\ / s,Ns ^ / l,s' 

,2 ^M 



V y s,Ns\ y i,s' 

+ Uis(w^G^) (W^G^) 0,/ 

'V J s,l\ J Ns,s' 



Summing over internal indices we find that 



Vi^sGgM^Gi,U^s',i 



=2«{0 



(e"W - 6o)2 (1 - e"(06(/))2 






-"(0 



-"(0 



-a{l) 



Ml) 



Ml) 



Ml) 



b{l 
b{l 
b{l 
b{l 
b{l 
b{l 



Gi{l) 
Gi{l) 
Gi{l) 
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All other terms in (|149|) vanish as Ng -^ co. Substituting ([L50D-(p!57[) into (|149|) gives the final 
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result: 
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